Entanglement entropy for the n-sphere 

H. Casini* and M. Huerta^ 
Centro Atomico Bariloche, 8400-S.C. de Bariloche, Rio Negro, Argentina 



Abstract 

We calculate the entanglement entropy for a sphere and a massless scalar field in 
any dimensions. The reduced density matrix is expressed in terms of the infinitesimal 
generator of conformal transformations keeping the sphere fixed. The problem is mapped 
^NJ ■ to the one of a thermal gas in a hyperbolic space and solved by the heat kernel approach. 

The coefficients of the logarithmic term in the entropy for 2 and 4 spacetime dimensions 
are in accordance with previous numerical and analytical results. In particular, the four 
dimensional result, together with the one reported by Solodukhin, gives support to the 
Ryu-Takayanagi holographic ansatz. We also find that there is no logarithmic contribution 
to the entropy for odd spacetime dimensions. 
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oT 1 Introduction 

£/ 

The entanglement entropy has proved to be of relevance within many research fields, in par- 
ticular, within the black hole physics [1] and quantum field theory contexts [2j [3] . It can be 
interpreted as a measure of the quantum correlations between two subsystems separated by a 
surface. By general considerations, the entanglement entropy Sy associated to a spatial region 
V in quantum field theory can be expanded as 

O : S v = ^-i£" (d_1) + - + die' 1 + go log[e] + S , (1) 

O 

where d is the dimension of the space, e is a short distance cutoff, the coefficients g>j are local 

functions of the edge surface dV and So is the finite part. This is due to the short range 
character of the ultraviolet quantum correlations, which makes the divergent part of Sy to be 
determined by the geometry of dV. In (pQ), the coefficients gi with i > are non universal 
in the sense that they depend on the regularization prescription. The non universal terms 
diverge with inverse powers of e. On the other hand, the logarithmic term with coefficient go 
is universal and then it has particular relevance since it is strictly related with the concerning 
theory. 

Here, we study this term for the entanglement entropy associated to a sphere of radius R 
in d spatial dimensions {{d + 1) dimensional Minkowski spacetime) and a massless scalar field. 
For spherical sets and when the considered field theory is conformal invariant, the entropy 
writes Sy = go\og(e/R)+non universal terms, and the coefficient go gives the only universal 
information contained in Sy\j. This coefficient has been related to the c-theorem in the context 
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1 Contributions to go depending on the radius appear for a sphere in non conformal theories [3]. 
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of the Maldacena duality [5]. In two spacetime dimensions, this coefficient is —(1/3) times the 
Virasoro central charge of the respective conformal field theory. For a scalar in two dimensions, 
it is go = — |. This is a very well known result which has been obtained using different 
approaches, and it is among the first exact results available in the field [6J. 

Fursaev and Solodukhin [7] have developed a generalized heat kernel expansion which is 
applicable to manifolds with general metric and conical singularities. This was used to evaluate 
divergent terms in the entropy in four dimensions using the replica trick. In [7J |S], it was used to 
obtain the value of go for black holes, which have horizon surfaces with zero extrinsic curvature. 
However, for flat space the contribution of the extrinsic curvature of dV is crucial in the 
geometric entropy. This contribution in four dimensions was obtained recently by Solodukhin 
in [9] . His approach consists in demanding conformal invariance of the most general geometric 
expression for the logarithmic contribution, and then using the Ryu-Takayanagi holographic 
ansatz [TO] (which is applicable to the special case of the entanglement entropy of SU(N) 
superconformal gauge theories) in order to calibrate a free coefficient. In particular, it gives 
go = 1/90 for the sphere and a scalar field. Recently, the same coefficient was calculated 
numerically in [TTj improving the numerical method of Srednicki [12] . The numerical result is 
also go = 1/90. An earlier paper questioned the validity of the application of the holographic 
ansatz for regions with boundary having nontrivial extrinsic curvature [13] . These doubts seem 
to be less justified in the light of the agreement of the results obtained by different methods on 
the value of go for the sphere. 

In this work we present the exact logarithmic coefficient for any dimensions computed 
analytically from first principles. In particular, the result for two and four dimensions agree 
with the results mentioned above. We also give the corresponding coefficients for the Renyi 
entropies. 

The entanglement entropy is derived from the vacuum state reduced density matrix py = 
£r_y|0)(0| on the set V by S = — tr (p v log py). In general, a density matrix, as a Hermitian 
and positive operator, can be written as the exponential of a Hermitian operator. We write 

pv = c e~ H . (2) 

This H is known in the literature as the modular Hamiltonian and it plays a central role in the 
algebraic approach to axiomatic quantum field theory [TJ]. There are only few cases in which 
% is exactly known. When V corresponds to half of a spatial plane, the modular Hamiltonian 
is proportional to the boost operator which keeps the wedge invariant [15]. This is a different 
way to express the result by Unruh: Taking the flow generated by the boost operator as a time 
flow, which coincides with the time meassured by constant accelerated observers, the vacuum 
state 02} in the region accessible to these observers is seen as a thermal one [IE] . 

For conformal theories and spherical sets the modular Hamiltonian is known to be pro- 
portional to the generator of conformal transformations which keeps the sphere fixed |17[ [T8] . 
This follows from the result for the half-space by a conformal transformation. Thus, in this 
particular case, H has the remarkable property of being local. 

We note that the relation ([2]) always admits a thermal interpretation of py. Consider an 
internal time r and the unitary evolution given by the operator py inside the algebra of the 
fields in V. Then, py becomes thermal with respect to the time translations operator 7i, with 
temperature rescaled to 1. Our calculation of the entropy in this paper is based on this idea. 
We found that for any dimension, the state of the massless scalar field inside the sphere can 
be mapped to the one of a thermal gas of free bosons at temperature 1 living in a hyperbolic 
space. Then, we calculate the coefficient go of the logarithmic term in ([[]) using the heat kernel 



approach in the evaluation of the corresponding partition function. 



2 The thermal gas 

For a conformal field theory and spherical sets, the modular hamiltonian is proportional to 
the generator K of infinitesimal conformal transformations which keeps the sphere invariant 

[mug 

% = 2nK . (3) 

The normalization of K is fixed to recover the Unruh thermal effect p ~ e~ 2irK , with K the 
boost operator, near the sphere boundary. 

In particular, for a massless real scalar field 4>(x), % is given by 



n = 2n I dV {R2 4R r2) U\x) + (V<f>(x) • V0(x)) - ^W 



(4) 



where the integration J dV = J drr^ d ~^ J dVt is taken over the volume of a (i-ball of radius R, 
and the integrand is the density j° of the Noether current corresponding to the one param- 
eter group of conformal transformations leaving the sphere fixed. We rewrite the quadratic 
expression (j4j) as 



n 



dV x / dV y (<k(x) N(x,y) ir(y) + <j){x) M(x,y) </>(y)) 



(5) 



in terms of the kernels M(x,y) and N(x,y). The tt(x) in (|3j) is the momentum canonically 
conjugated to the field <ft(x), which satisfy the standard equal-time commutation relation 



[n{x),<f>{y)] = -iS(x-y) . 
From (j4j) we identify the explicit expressions for M and N, 



with f(r) 



M 
N 

(R 2 -r 2 ) 



-2n5 d (x - y)[f(r)A + V/(r) • V + ^J(A/(r))] 



(6) 



(7) 



2n6 d (x-y)f(r), (8) 

and A = Tz-f§-(r d ~ l jl-) + -^Aq the Laplacian in d dimensions. After 



2R "" iivl " ~~ r 1 *- 1 d r V d r 

a Bogoliubov transformation of the fields and momentum variables, the Hamiltonian takes a 
diagonal form and the density matrix ()2]) writes (2J [20] 



Pv 



na 



(9) 



in terms of independent normalized boson operators a n , a^ and the energy eigenvalues e n . Eq. 
(jUJ) corresponds to a thermal state of independent bosons with energy e n and temperature 
T — 1. We note that the fact that the considered theory is free is what makes possible the 
mapping to a thermal gas of independent scalar particles. 

The square of the energies, e^, is given by the eigenvalues of the product kernel NM [2] 



Qe 2 Q~ 1 = NM = -47i 2 5 d {x - y) 



fA + ffd r + ^2^/(A/) 



(10) 



for some invertible operator Q. Taking into account OH]) and fTTTJj) is then possible to map the 
original problem of calculating the entanglement entropy for a scalar free field reduced to a 
spherical set, into a new one consisting on the evaluation of the entropy for a thermal gas of 
scalar particles with temperature T = 1 and Lagrangian 

C=^(-df-NM)(f>. (11) 

Indeed, the field equation for this theory gives the same spectrum as (ITUj) for the energies. The 
relation between the entropy and the corresponding partition function is the standard one for 
a thermal state 

S= [1-Pdp] log Z(p)\p=i ■ (12) 

In the next section we calculate the effective action W = — log Z(j3) in terms of the heat kernel 
associated to the NM differential operator. 

3 Heat kernel approach 

Consider a general second order differential operator D of the Laplace-Beltrami type 

D = -(g^d fl d u + a°d a + b), (13) 

where g^ v can be interpreted as the inverse of a metric tensor. Following standard procedures 
], D can be rewritten in the form 



D = -(g^V^ u + E), (14) 

where the covariant derivative V = V^ + w contains both Riemann V'^ and gauge w parts. 
We may express E in terms of a M and b in ffT3|) 

E = b-g^{d^w u + w u w p -w a T^) 1 (15) 

where the connection w^ and its Christoffel symbol T are given by 

ws = iflw(o" + ^T^), (16) 

K» = 9 up \{d,g vp + d v g^-d p g^). (17) 

The expression (Il4"|) of the differential operator is more suitable for the application of the heat 
kernel in the evaluation of the partition function (TT2l) . 

According to (ITUj) and (TT3"|) . the operator D = NM gives in spherical coordinates 

9^ = a 2 f 2 dia g( 1 ' r<2 ' r2sin6l i » r 2 sin 6>i sin 6*2 , ... , r 2 sm6l ... sin^_ 2 ) , (18) 



4tt 2 / 2 
(d- l)r 



0, ... , , (19) 



r 2 -R 2 '■ 

E (d-l) 2 7r 2 . (20) 

The field strength d^vjy — d v w^ of the connection w p is zero for all dimensions. Therefore, w p 
can be removed by a "gauge" transformation which leaves invariant the eigenvalues, and makes 
hermitian the differential operator NM 

Wfl ^w'^ = w^- d^r) = , 0^0' = e h(r) , (21) 
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with h(r) = &=& \og(R 2 - r 2 ). 

The sectional curvature of the metric (fl8|) is constant and negative 



*5^,«w __, 2 = _ te2 _ 



U a U a V/3vP — (u a V a ) 2 



for any vectors u and t>. Here -R^ pCT is the curvature tensor. Thus, the metric corresponds to 
a d-dimensional hyperbolic space H^. Note that the sectional curvature is independent of the 
dimension. 

Thus, we can consider just the operator 

D = -^V[f ] V^ 1 - (d - 1) V (23) 

in this hyperbolic space. It is notable that the constant E = (d — l) 2 ir 2 exactly eliminates the 
gap in the spectrum of the laplacian in Hd, leaving D gapless. 

Then, we consider the partition function in the thermal theory corresponding to (TTTT) . In the 
Euclidean approach, and once the time coordinate is compactified (in our case we take (3 = 1 
at the end of the calculation), the (d + 1) space-time is S\ x Hd, where S\ is a circle of size (3 
and Hd is a d- dimensional hyperbolic space with sectional curvature — k 2 . 

The partition function Z(f3) is related to the heat kernel [21] by 

1 r°° dt 
lOgm = 2j J K{t) - (24) 

The heat kernel K(t) for a differential operator O is defined as 

K(t) = Tr(e- to ) = f d d Xy/gK(t,x,x) , (25) 

with K(t,x,y) = (x\e~ ot \y) . For the operator — d 2 + D on the product manifold S\ x Hd we 
have 

K(t) = K Sl (t) x K 3 (t) . (26) 

Here, 



DC 



V47rt ^— f 

v 71=1 

is the regularized heat kernel of — d 2 for the circle, in which the f3 — > cxd (zero temperature) 
limit has been subtracted. We also have 

K 6 (t) = K Hd (t)e E \ (28) 

where Kn d is the heat kernel of the hyperbolic space defined by a pure Laplace-Beltrami 
operator (E = in flHJ)), and E is the constant given in (|2"U|) . 

Since the hyperbolic space is homogeneous the volume factorizes in the heat kernel, 



K Hd {t) = K Hd (t,x,x) ld d x^~g. 



(29) 



In order to regularize the entropy, we have to cut off the ultraviolet modes near the surface of 
the sphere. This is done integrating up to a radius r < R — e. Remarkably, this is seen as an 
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Table 1: Coefficients of the logarithmic term in the entanglement entropy associated to a sphere 
for even spacetime dimensions. 



infrared regulator in hyperbolic space, since r — >• R corresponds to an infinite geodesic distance 
limit. Expanding yfg in powers of e = R — r, the volume integral can be written as 



d d x y/g = I dQ dr 
The coefficient 



„d-l 



(J dQ) f (d ) R d (d) 



(2nf(r)Y (2nY I &-17ZT + - + « ^(e/R) + ... 



d-l 






EK 



d+j_ 



(d + j-iy. 



j=0 



(2j)\\{d-3-l)\j\ 



(30) 



(31) 



is a multiplicative factor in the coefficient of the logarithmic term of the entropy. From the 
Taylor expansion in ( 130]) it follows that there is no logarithmic contribution to the entropy in 
odd spacetime dimensions since (quite mysteriously) q$ 



(d=2m) 



in (131)1 . Curiously, this is the 



only coefficient which vanishes in the series. 

On the other hand, for odd spatial dimensions, d = 1m + 1, the heat kernel Kn d is [22) 



K Hd {x,y,t) 



1 



d 



I 1 \m L,2m+1 

(2vr) m (47rA; 2 t)( 1 /2) \ksmhkpdp 



e -k 2 mH-^ 



d = 2m + l 



(32) 



with p the geodesic distance between x and y, and — k 2 = —4n 2 the sectional curvature. From 
this expression, it follows that for m > 0, K Hd (t, x, x) takes the general form 



KHd{t ^ x) = E^L e ^-^ 



(Ant) 



(33) 



with Pd(x) a polynomial. The exponential is due to the gap in the spectrum. We have Pi(t) = 1, 
and we write 

(d-3)/2 

p d (x)= J2 a f xi > ( 34 ) 

i=o 



for d > 3. The coefficient aS are rational numbers. Then, from (128)1 the exponential factor 
cancel, and we have 



(47rt)2 y 



Finally, the entanglement entropy is obtained combining f[T2~)) . ([24)1 . (13TJ)) and 
coefficient of the logarithmic term is go — — 1/3 for d — 1, and 



(35) 
The 



9o 



(d) n (d-3)/2 



2^r[d/2] ^j ' 



(d+l)/2 



(d) 

g 

2^(d-2)H 



Em 



>m Jrf) 



d+l-2i 



/=2 



2 

2 2l U 
(27) 



■B 



2/ 



(36) 



for d odd greater than one. Here the coefficients q$ and aS are given by (13TT) and ([31 
respectively, r[x] and ([x] are the standard Gamma and Riemann zeta functions, and B n are 
the Bernoulli numbers. The results for the coefficients for spacetime dimensions up to d+1 = 14 
are given in table 1. 

Another class of interesting information measures are the Renyi entropies 

S„ = i^l. (37) 

1 — 71 

These are usually computed using the replica method, and then the limit lining! S n = S is 
taken to obtain the entropy. Here we can express S n in terms of the partition function 

S n = ~^— (log(Z(n/3)) - n log(Z(/3))) . (38) 

1—71 

Then the logarithmic coefficient in the Renyi entropies (S n = golog(e/i?)-|- non-universal 
terms) is 

jn) _ % sr r-iv+^w [i-iy.ti2i W -1 . , 

i/o ~~ d+i , , . .. / j \ 1 a d+i~2i / on . / 1 N 27 _i • w»; 

2— (d-2)!! ^ ^^ ( 2 ! (n-l)n 2 ' 1 

These are rational functions of n which converge to non zero values for n — > oo. More explicitly, 
for d — 1, 3 and 5 we have 

<T = ~ , d=l, (40) 

w = (n-t-1) (»' + !) , = 3 

y ° 360 n 3 ' V ; 

= _(„ + l)( 3n » + l)(3 n ' + 2) 

y ° 30240 n 5 ' v ; 

4 Final Remarks 

Solodukhin [9] provided a general formula for the logarithmic coefficients in the entropy for 
d + 1 = 4 dimensions, in terms of conformal anomaly coefficients. In that work there is a free 
coefficient which was calibrated using the Ryu-Takayanagi holographic ansatz for the entan- 
glement entropy of super conformal gauge theories. Our result go = 1/90 for the logarithmic 
coefficient in four dimensions agrees with the Solodukhin formula (and also with the numerical 
results [II]). We can thus turn the argument around, and conclude that our analytical result 
gives support to the holographic ansatz. 

The mapping of the problem on a sphere in flat space to one in a hyperbolic space is unex- 
pected. However, it is known that the near horizon geometry of several black hole spacetimes is 
AdS space (which is spatially hyperbolic). Emparan also discusses the possibility to define the 
entropy in Rindler space approaching from a spatially hyperbolic space [23]. A map relating 
the entropy in S* 1 x H & to a sphere was also used in the context of the AdS-CFT duality in [5] . 
It would be interesting to explore if the spectral problem of reduced density matrix for a sphere 
and different conformal theories, in particular the free massless fermion, can also be mapped 
to equivalent problems in Hj. 

Finally, we remark a curiosity of the result for two and four dimensions. In these cases the 
entanglement entropy of the sphere is equivalent to the one of a gas with the flat space entropy 



density, but living in hyperbolic space (the volume has to be calculated taking into account 
the curved geometry). In higher dimensions however, there are curvature corrections to the 
entropy density. 

Note added 

After the circulation of the present paper as a manuscript Dowker [2Ij presented a calculation 
of the coefficient go for a sphere and a scalar field in any dimensions by a different method. His 
results agree with the ones in this paper. Also, Solodukhin [25] has shown the same logarithmic 
coefficients arise for the extremal black hole. 

Acknowledgments 

This work was partially supported by CONICET and Universidad Nacional de Cuyo, Argentina. 

References 

[I] L. Bombelli, R. K. Koul, J. Lee, and R. D. Sorkin, Phys. Rev. D 34, 373383 (1 986); 
C. G. . Callan and F. Wilczek, Phys. Lett. B 333, 55 (1994) |arXiv:hep-th/9401072] . 

H. Casini and M. Huerta, J. Phys. A 42, 504007 (2009) |arXiv:0905.2562l [hep-th]]. 

P. Calabrese and J. L. Cardy, J. Stat. Mech. 0406, P002 (2004) [arXiv:hep-th/0405152] ; 
P. Calabrese and J. Cardy, J. Phys. A 42, 504005 (2009) [ arXiv:0905.40l3l [cond-mat.stat- 
mech]]; J. L. Cardy, O. A. Castro-Alvaredo and B. Doyon. larXiv :0706.3384 [hep-th]. 

M. P. Hertzberg and F. Wilczek, larXiv: 1007.0993 [hep-th]. 

R. C. Myers and A. Sinha, Phys. Rev. D 82, 046006 (2010) jarXiv: 1006. 12631 [hep-th]]. 

See for example, C. Holzhey, F. Larsen and F. Wilczek, Nucl. Phys. B 424, 443 (1994) 
|arXiv:hep-th/94031~0~8J . 



[2] 
[3] 

[4] 
[5] 

[6] 

[7] 



[9 
[10 

[11 
[12 



S. N. Solodukhin, Phys. Rev. D 51, 609 (1995) |arXiv:hep-th/9407001|; D. V. Fursaev and 
S. N. Solodukhin, Phys. Lett. B 365, 51 (1996) |ar Xiv:hep-th/9412 020|; D. V. Fursaev and 
S. N. Solodukhin, Phys. Rev. D 52, 2133 (1995) | arXiv:hep-th/9501127] . 

R. B. Mann and S. N. Solodukhin, Nucl. Phys. B 523, 293 (1998) |arXiv:hep-th/9709064| ; 
S. N. Solodukhin, Phys. Rev. D 57, 2410 (1998) | arXiv:hep-th/9701106] . 

S. N. Solodukhin, Phys. Lett. B 665, 305 (2008) [arXiv:0802.3TT7l [hep-th]]. 



S. Ryu and T. Takayanagi, JHEP 0608, 045 (2006) |arXiv:hep-th /0605073|; T. Nishioka, 
S. Ryu and T. Takayanagi, J. Phys. A 42, 504008 (2009) [arXiv: 0905 .09321 [hep-t h]] ; S. Ryu 
and T. Takayanagi, Phys. Rev. Lett. 96, 181602 (2006) |arXiv:hep-th /0603001~] . 

R. Lohmayer, H. Neuberger, A. Schwimmer and S. Theisen, Phys. Lett. B 685, 222 (2010) 
|arXiv:0911.4 283 [hep-lat]]. 



M. Srednicki, Phys. Rev. Lett. 71, 666 (1993) |arXiv:hep-th/9303048|. 



[13] A. Schwimmer and S. Theisen, Nucl. Phys. B 801, 1 (2008) |arXiv:0802. 10171 [hep-th]]. 

[14] R. Haag, Local Quantum Physics: Fields, particles, algebras, Springer (1992); H. Borchers, 
J. Math. Phys. 41, 3604 (2000). 

[15] J. J. Bisognano and E. H. Wichmann, J. Math. Phys. 17, 303 (1976). 

[16] W. G. Unruh, Phys. Rev. D 14, 870 (1976). 

[17] P. D. Hislop and R. Longo, Commun. Math. Phys. 84, 71 (1982); 

[18] H. Casini, Class. Quant. Grav. 25, 205021 (2008) jarXiv:0804.2~T82l [hep-th]]. 

[19] H. Casini and M. Huerta, Class. Quant. Grav. 26, 185005 (2009) [arXiv:0903.5284l [hep- 
th]]. 

[20] I. Peschel, J. Phys. A: Math. Gen. 36, L205 (2003) |ajXiv:cond-mat/021263i] . 

[21] See for example: D. V. Vassilevich, Phys. Rept. 388, 279 (2003) [arXiv:hep-th7 0306138|. 

[22] A. Grigor'yan, M. Noguchi, The heat kernel on hyperbolic space, Bull. Lond. Math. 
Soc. 30, (1998) 643. A. Grigoryan, Upper bounds on a complete non compact manifold, J. 
Funct. Anal. 127, (1995) 363. A. Debiard, B. Gaveau, E. Mazet, Theoreme de comparison 
in geometrie riemannienne, Publ. Kyoto Univ. 12 (1976) 391. 



[23] R. Emparan, Phys. Rev. D 51, 5716 (1995) |arXiv:hep-th/9407064| 
[24] J. S. Dowker. la7Xiv:1007.3865 l [hep-th]. 
[25] S. N. Solodukhin. TarXiv: 1008.4311 [hep-th]. 



